Fokker-Planck equations and density of states in disordered quantum wires 
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We propose a general scheme to construct scaling equations for the density of states in disordered 
quantum wires for all ten pure Cartan symmetry classes. The anomalous behavior of the density 
of states near the Fermi level s = for the three chiral and four Bogoliubov-de Gennes universality 
classes is analysed in detail by means of a mapping to a scaling equation for the reflection from a 
quantum wire in the presence of an imaginary potential. 

PACS numbers: 71.55.Jv, 71.23.-k, 72.15.Rn, 73.23.-b 



S3 
S3 

C/3 



I 

S3 
O 

o 



> 



O 
O 



T3 

S3 
O 

o 



X 



I. INTRODUCTION 

Statistical properties of energy levels and wavefunc- 
tions in disordered electron systems are believed to be 
determined by, first of all, the fundamental symmetries 
of the Hamiltonian. Such a connection is best estab- 
lished in the case of "zero-dimension", i.e., for finite 
size systems with a large dimensionless conductance g, 
where a description in terms of random matrix theory 
is valid. Usina-the link ta-Cartan's classification of sym- 
metric spacesjd Zirnbauero and Casellea have pointed out 
that there exist only ten possible random matrix theories, 
whose form follows directly from the geometrical charac- 
teristics ( "roots" ) of the corresponding symmetric space 
in Cartan's table. These ten random matrix theories am 
divided into three standard classesja three chiral classes,El 
and four Bogoliubov-de Gennes (BdG) classesjj the sub- 
division in each class depending on the presence or ab- 
sence of time-rewersal symmetry (TR) and spin-rotation 
invariance (SR)j3 see Table |. 

The rational for the Cartan classification is believed to 
transcend "zero-dimension" and has been applied to the 
construction of effective theories (e.g., non-linear-sigma- 
models) for higher dimensional disordered systems. The 
standard classes are thus believed to be appropriate to 
the problem of an electron maving in a random potential, 
without further symmetries,!! irrespective of dimension- 
ality. The chiral classes are appropriate to the case when 
the disorder is purely "off-diagonal" p as is the case, e.g., 
for the lattice .tandom flux model in quasi-QaCp-and two- 
dimensions ,t!30 the random hopping model and for 
random XY spin chains .Eil The BdG classes refer to sys- 



tems with superconducting correlations ,c£l and— were ar- 
gued to be valid for vortices in superconductors, E30 dirty 
unconventional superconductorsjL3o and (in the case of 
broken time-reversal symmetry) normal metals in prox- 
imity to a superconductor. The BdG classes have also 
been argued to be of some relevance to problems in statis- 
tical mechanics such as random bond Ising and network 
modelsM 

The geometric structure of the symmetric spaces not 
only determines the level statistics in zero-dimension, it 



also determines the form of scaling equations for the 
transmission eigenvalues Tj and reflection eigenvalues 
Rj = 1 — Tj (j = 1, . . . , N) in a wire geometry, where a 
disordered sample of Vpigth L is connected to two ideal 
leads (see Fig. [^a).dEjO Here the reflection eigenvalues 
Rt, . . . , Rn are the d-fold degenerate eigenvalues of the 
matrix product r'r, where r is the iVd-dimensional ma- 
trix of reflection amplitudes for back-reflection into the 
same lead (see Fig. [j]a) and Nd is the number of prop- 
agating channels at the Fermi level. The degeneracies d 
are listed in the Table | for the 10 symmetry classes in 
Cartan's classification. Parameterizing 



Ri 



tanh Xn 



the scaling eq 
takes the for" 1 ' 1 ~' 



J ~ 3 ' 

n for the distribution P{x\, 

N 



dL 2j£^dxi dxi 



(1) 
,x N ;L) 

(2) 



where I is the mean free path and L is the length of 
the quantum wire. For the standard and BdG symmetry 
classes, the Xj are positive random variables, < Xj < 
oo, and the Jacobian J and the numerical constant 7 
read 

N 

j=yi smhm! ( 2x j) n n sinhm ° ± ( 3a ) 

j = m (N -l) + mi + l, (3b) 

while for the chiral classes the Xj can take values on the 
entire real axis, -co < Xj < 00, and one has 



j=nn sinhmo (^- 

j=i k<j 

7=| [m (N - 1) + 2] . 



Xk), 



(4a) 
(4b) 



Here m; and m are the multiplicities of the so-called 
long and e*dinary roots for the corresponding symmet- 
ric space,c3 respectively, see Table || (For the chiral 
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FIG. 1. A quantum wire (hashed marked region) of length 
L connected to ideal leads. The reflection matrix at energy e, 
r(e), is the matrix of reflection amplitudes for reflection into 
the same lead, (a) A two lead geometry is used to compute 
the conductance, (b) A one lead geometry is used to compute 
the DOS whereby the wire is closed by a perfectly reflecting 
wall on the right. 



classes we only consider the case when the nonuniver- 
sality parameter 77 equals unity, see Ref. |2(|) Fpx the 
standard classes, Eq. (||) waSrderived by DorokhovjHZl and 
Mello, Pereyra, and Kumar;c3 derivations for the chiral 
and BdG symmetry classes can be found in Refs. ^9] and 

E- 

In this note, we show how the same structure also de- 
termines scaling equations for the density of states (DOS) 
u(e) in a quantum wire of infinite length (L — > 00) at en- 
ergy e. Our work builds on previous work for the chiral 
classes, where such .a. scaling equation was derived using 
a different method.Ej In addition, we present an exact 
solution for the DOS, something that could not be done 
in Ref. |lj 

While v(e) is a non-singular function of energy for the 
three standard classes, singular behavior is expected near 
e = for the remaining seven symmetry classes, as e = 
is a special point there. (The energy e = corresponds 
to a point of particle-hole symmetry in the chiral and 
BdG symmetry classes; it is the band center for lat- 
tice models with random hopping, or the Fermi energy 
in the case of the BdG classes.) Indeed, for the chiral 
classes, v{e) was found to depend sensitively on the par- 
ity of the channel number N for e close to zero, showing 
the v(e) oc l/|eln 3 e| divergence characteristic of pure 
one^dimensional systems with chiral symmetry for odd 
AT,H while v{e) oc le" 1 '- 1 lne| for even NH For the BdG 
classes in the presence of spin-rotation invariance (classes 
C and CI; we refer to the symmetry classes by Cartan's 
symbol for the symmetric space corresponding to their 
Hamiltoniaas.. see Table ||), a suppression u(e) oc e" 1 '^ 1 
is expectedJiSH In both cases, the characteristic energy 
scale for the DOS singularity is ~ Hvf/N 2 £, vf being 
the Fermi velocity. This distinguishes the singularity 
in v(e) for quantum wires from its counterpart in zero- 
dimension, where the characteristic energy scale is the 
mean level spacing, which goes to zero as the system size 
is increased. No exact results are known for the multi- 
channel quantum wires from BdG symmetry classes in 
the absence of spin-rotation invariance (classes D and 
Dili; see however Refs. |34|-|36| for asymptotic results in 
one- and two-dimensions). The general scheme that we 
present here fills this gap, and premdes a unified frame- 
work for all ten symmetry classes. Ej 
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TABLE I. Classification of symmetry classes according to 
Cartan's table. The symmetry classes are defined in terms 
of the presence or absence of time-reversal symmetry (TR) 
and spin-rotation invariance (SR), and in terms of the other 
fundamental symmetries of the system (standard, chiral, or 
BdG). The table further contains the multiplicities of the or- 
dinary and long roots m and m; audi the degeneracy d of the 
reflection/transmission eigenvaluesEfj In the seventh column 
we denote, using Cartan's notation (chapter X of Ref. [I]), the 
symmetric spaces associated to the transfer matrix group M 
of the quantum wire; these are different from the symmetric 
spaces associated to the microscopic Hamiltonian Ti of the 
quantum wire, which are listed in the last column. Following 
the convention of Ref. [| we refer to the 10 symmetry classes 
by Cartan's symbol for the symmetric space of their Hamil- 
tonians. For standard and chiral quantum wires without TR 
the degeneracy d of the reflection eigenvalues is 2 (1) with 
(without) SR. 
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II. MAPPING TO QUANTUM WIRES WITH 
ABSORPTION 

Our construction of a scaling equation describing the 
DOS makes use of a mapping to a quantum wire with 
absorption. Absorption is described by the addition of 
a spatially uniform imaginary potential iu>, lu > 0, to 
the Hamiltonian. In the presence of absorption, a scaling 
equation for the reflection eigenvalues Rj, j = 1, . . . ,N, 
of r(ibj)^ r(ioj) can be derived in the same way. as the scal- 
ing equation (j|) in the absence of absorption,E3c3 noting 
that upon addition of a thin slice of length SL to a dis- 
ordered wire of length L, the reflection eigenvalues Rj 
change by a small amount SRj. The change SRj comes 
from two statistically independent contributions: (1) dis- 
order scattering in the added slice of width 5L, and (2) 
absorption in the same slice. These contributions are 
statistically independent, since to linear order in 5L one 
can neglect processes that involve both disorder scatter- 
ing and absorption. The second process gives rise to a 
change 



R i 



-4oj<5L j vp 



, N, 



(5) 



independent of the disorder. The first process depends on 
the disorder type and symmetries, and is already present 
in the scaling equation (^). Combination of the two pro- 
cesses then results in the replacement of Eq. (|J) by 



dPiu, 
dL 



N 



^ dxj \v F 2j£ dxj J 



3=1 



(6) 



where we used the parametrisation ([l]) for the Rj, and 
Piu) {x\ , . . . , xjv ; L) represents the joint distribution of the 
reflection eigenvalue in the presence of the imaginary po- 
tential iijj. Equation (JsJ) was derived in Refs. |3^j3^ for 
the standard symmetry classes. The initial condition 
at L = for Piu(x\, . . . , xn; L) still needs to be speci- 
fied. For our purposes, it is advantageous to set Rj = 1, 
j = 1, . . . , N, at L = 0, corresponding to a wire that is 
closed on one end, see Fig. 

The scaling equation (|6|) has the stationary (indepen- 
dent of L) solution, 



JY 



Piuj;st (^l ) • 



v 1 j=i 



-a cosh 2xj 



(7) 



where a = j£lu/vf and the normalization constant Z(a) 
ensures that the stationary solution Piuj-st is normalized 
to one. This solution represents the joint probability 
density of the reflection coefficients for the absorbing 
medium. We will demonstrate below that the DOS in 
an infinitely long disordered wire without absorption can 
be computed as a logarithmic derivative of the normaliza- 
tion constant or "partition function" Z{a) with respect 
to the dimensionlcss imaginary energy a. 



The main idea of our method is to make the ana- 
lytical continuation iuj — » e. With this analytical con- 
tinuation, the single-energy product r(iujYr{iuj) of re- 
flection matrices in the presence of absorption maps to 
the two-energy product r(— £yr(e) without absorption.c2l 
The latter product is a unitary matrix since the reflection 
matrix in the geometry of Fig. |l|b is unitary by construc- 
tion. The TV independent eigenvalues of r{—eyr{e) can 
thus be written as exp(2i(f)j), j = 1, . . . ,N. When the 
length L of the disordered wire is increased, the <f)j 's per- 
form a Brownian motion, governed by a Fokker-Planck 
equation that is the analytical continuation of Eq. (^), 



dP e _ 9 ( 2e , 2 „ ; , 2 



where 



(8b) 



k<j 



for the standard and BdG classes, and 



/v 



II II 



(8c) 



k<j 



for the chiral classes, where 7 is defined in Eqs. (ph and 
([I|), respectively. This equation was derived by different 
methods in Ref. |3l] for the chiral classes and in Ref. ^l] 
for the standard classes. The scaling equation (||) can be 
considered either on the full real axis —00 < <j)j < 00, 
or on the interval — 7r < <fij < tt. In the latter case 
Eq. (||) has a well-defined stationary solution P e . st ({(j)j}) 
describing the distribution P e ({4>j}; L) in the limit L 
00. However, this solution is not of the simple form (Q).t3 
In fact, no representation of P e;s t in terms of elementary 
functions could be found other than for the caseJV = I, 
where explicit integration of Eq. (^) is possibles 

The DOS is calculated via the "node-counting theo- 
rem", which states that, for the half-open quantum wire 
of Fig. 0b, the mean number of states Af(e) in the energy 
interval (— e, e) and per unit length is proportional to 
the sum of the phase derivatives d<j>j / '<9L.E-ac3 For tech- 
nical reasons, we replace d(j)j/dL by the derivative of the 
counting function — Im lnsin^-HO), which has the same 
average slope as (f>j. We thus find 



me) 



de'v{e') 



d" d 



= ~~ W < Im lnsin (^ + M ))eJ 



7T ^ 8L 



(9) 



Note that the positive infinitesimal iO in the argument of 
s'm(c/>j + iO), which is needed to select the correct branch 
of the logarithm, is compatible with the analytical con- 
tinuation iuj — ► e used to obtain the scaling equation m). 
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In Eq. (||), the brackets (. . .) s> l denote an average with 
the probability density P e that solves Eq. (|§|) on the full 
real axis — oo<cf>j < oo at length L of the quantum wire. 
Integrating by parts, with the help of Eq. dq), we find 



• A/ '( £ ) = --E *h 

3=l J 



Im [cot(0j + iO)] 



d 



2de 

TTVf 



N 



5^(Im COt(0j + i0)) e 



(10) 



3=1 



Now, in the thermodynamic limit L — > oo, the disorder 
average (. . .) E can be performed with the stationary dis- 
tribution P £ - S t for the interval — n < 0j < 7T. Hence the 
DOS per unit length is given by 

. . 1 dJ\f d d . , , . , . 

= - -^ = - — £ ^ ( I- cotfo + ,0)) £ , (11) 

i=i 

where we have used the relation 



v{e) = v{-e) 



(12) 



that holds for each disorder configuration in the chi- 
ral and BdG classes. In the standard classes, Eq. Jl2] ) 
only holds in the thermodynamic limit and for suffi- 
ciently small e (smaller than the energy scale vp /£, where 
nonuniversal effects appear and the validity of the scaling 
equations breaks down). 

Instead of performing the disorder average with the 
stationary solution P e . st of Eq. (||), we rely on the solu- 
tion Piuj-st for the stationary distribution of the reflection 
eigenvalues in the presence of the imaginary potential ito, 
see Eq. (Q) . From the relation tanh 2 Xj = e 2z ^ j we first 
observe that 



Im (cot(0j +i0)) £ 



Re (cosh 2xj 



(13) 



where (. . .) iL0 stands for averaging with the stationary 
distribution Pj w;s t- Second, we note that the average 
— ^2 m {cos\i2x m ) is given by the logarithmic derivative 
of the partition function Z(a), thus obtaining the DOS 



f \ d T> 9 

v(£) = Re — 

ttvf oa 



a— In Z(a) 
oa 



(14) 



This relation between the DOS and the partition function 
Z[a) for each of the Cartan symmetry classes is the key 
result of this paper. 



III. CALCULATION OF THE DENSITY OF 
STATES 

As we have seen in the previous section, the computa- 
tion of the DOS in the thermodynamic limit reduces to 
the computation of the "partition function" 



Z(a)= / dxi...dx N \J({x j })\e~ a Z'?=i 



cosh 2X4 



(15) 



For the standard and BdG classes, the integration over 
the coordinates Xj is restricted to the positive real axis, 
< Xj < oo, j — 1, . . . , N, while for the chiral classes, 
the integration extends over the entire real axis. 

For all ten symmetry classes, a change of variables 
xj — > Xj can be made after which Z(a) can be written in 
the form 



N 

Z oc / d\x . . . dX N Y[ J| l A J _ A 

Jc 3=1 k<j 



(16) 



The relation between the Xj and the Xj , the weight func- 
tion w(X), and the integration interval C depend on the 
symmetry class and will be specified below. 

A general method for calculation of integrals of the 
type (|l(]) is described in chapters 5 and 6 of Mehta's book 
on random matrix theory.Q Here we mention the results 
from Ref. |], up to a proportionality constant that does 
not depend on a. For m = 2 the result of integration in 
Eq. (jig) is proportional to the determinant 



Z oc Det [/ mB ]" = i , 



fn 



dXw(X) p m (X)p n (X), 



(17a) 
(17b) 



where the polynomials {p m (A)} form a linear indepen- 
dent set of polynomials up to degree N — 1. For rn = 1 
and even N we haveu 



Z OC J Det [fmn] 



N 

™»lm.n=li 



(18) 



f mn = J^dXdX' e(X- X')w(X)w(X')p m (X)p n (X'), 

where the sign function e(A) equals 1 for positive A and 
— 1 for negative A. For odd N, the matrix / must be 
supplemented by an additional row and column 



Z oc 



\f Det [fmn}Zt=V 



dXw(X)p n (X). 



fn,N+l — — JiV+l,n — 

For m = 4, the result is 
z K \l Det 



fmn = J^dXw(X)[p m (X)p' n (X) -p n (X)p' m (X)]. 



(19a) 
(19b) 

(20a) 
(20b) 



where the prime stands for the derivative with respect to 
A and {p«(A)} form a linear independent set of polyno- 
mials up to degree 2 N — 1 . 
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A. Standard classes 



(26b) 



For the standard universality classes, the partition 
function Z(a) takes the form of the normalization into* 
gral of the Laguerre ensemble of random matrix theoryO 
if we set A = sinh 2 x. In this case, the integration range 
C = (0, oo) and the weight function w(X) is given by 
cxp (— 2aX — a). Integrating over the Xj, one finds that 
the partition function Z(a) has a particularly simple de- 
pendence on a, 



Z(a) oc e- aN a^ N ' 2 . 



(21) 



Using Eq. (jlj), one quickly verifies that this corresponds 
to a constant DOS, 



Nd 

v{e) = v = . 

irvp 



(22) 



(The DOS can acquire a non-universal energy depen- 
dence on the scales s much larger then the inverse scat- 
tering time vp/£, where the applicability of the scaling 
equation breaks down.) 

As a matter of fact, for large a, one obtains the re- 
sult (|Tj) for all ten symmetry classes, as one verifies by 
substitution Xj — > a -1 / 2 ?/., in Eq. ( [T5| ) and subsequent 
expansion for large a. This confirms that the DOS v{e) 
approaches the standard value vq for energies far away 
from zero for all symmetry classes. 



B. Chiral classes 

For the chiral classes, we arrive at the standard form 
@ taking A = exp (2x), C = (0, oo), and 



;(A) = A-^ 2 cxp [-^(A + A- 1 )]. 



(23) 



When TR is broken, say by a magnetic field (m = 
2; class AIII), we choose in Eq. ([17]) the polynomials 
p„(A) = A™" 1 so that 



Z(a) oc Det [Aj m _ n[ (a)]^ , 



(24) 



where Kj(a) is the modified Bessel function of integer 
order j. For small values of the dimensionless imagi- 
nary energy a one can substitute Kq(o) — » —In a and 



K m (a) 



l (m- 1)!, m = 1,2, 



to obtain the 



leading behavior of Z (to logarithmic accuracy). For odd 
./V the determinant in Eq. (E4) is expanded according to 

Z(a) oca- (jy2 - 1 )/ 2 [lna + 0(l)]. (25) 
This gives the asymptotic DOS 

v ^ = i rS \ I ' < er « 1, (26a) 
\£T In (er) 

where we introduced the time scale 



It is the time needed to diffuse through a wire of length 
~ N£, or, equivalently, the mean DOS in a segment of 
the wire of length Nl. (The length scale Nt characterizes 
the crossover between the regimes of diffusive dynamics 
and of localized or critical dynamics.) For even A the 
determinant in Eq. ( p4| ) is expanded according to 

Z{a) oc a"* 2 ' 2 [1 - \ N 2 a 2 In 2 a + 0{a 2 In a)] . (27) 
This gives rise to a different asymptotic DOS, 



u(e) = 7Wo|£Thi (et)\, < et 1. 



(28) 



Observe that the asymptotic results (26a) and (|2q) do 
not depend explicitly on N when the energy is measured 
in units of 1/t. The leading dependence on energy in 

§s. (26a) and ( p8[ ) has been derived previously in Ref. 
in an approximation that leaves the prefactor -kvq un- 
specified. 

When TR symmetry is present but SR symmetry is 
broken by spin-orbit coupling (m Q = 4, class CII), we 
choose in Eq. (]2(]) the polynomials p2«-i(A) = A 2W_ ™ + 
A"" 1 and p 2 n (A) = \ 2N ~ n - A n_1 . With such a choice 
the square root of the determinant of a 2N x 2N matrix 
in Eq. ( po|) reduces to the determinant of a N x N matrix, 
from which we obtain 



iJV 



Z(a) oc Det [f nm \ n , m=1 , 
f nm = (2N -n-m + l)K\ n _ m \(a) 
+ (n - m)K 2 N-n-m+i(a)- 



(29) 
(30) 



For odd N, the leading term in the asymptotic expansion 
of Z(a) for small a is of the form (^5|), hence the small-e 



asymptote of the DOS is given by Eq. (26a), just like in 
the case of broken time-reversal symmetry. For even N 
the determinant in Eq. ( p9| ) is expanded as 

Z(a) oc a~ N2 [l + cia 2 + c 2 a 4 In 2 a + C(a 4 In a)] , (31) 

where c 2 = jgiV 2 (iV + l) 2 and ci is another A^-dependent 
coefficient, which does not affect the asymptote of the 
DOS. We thus obtain 



K £ ) = ^o(l + ^) 2 |(er) 3 ln( £ r) 



< er < 1. 



(32) 



In the presence of both TR and SR symmetry (m Q = 
1, class BDI) and with the choice of the polynomials 



PnW = A' 



1, 



, A, we obtain for even A 



Z(a) oc y Det [f mn ]Z,n=v 



(33) 



OO OO 



f mn = / dx / dy(xy)-3^ +L >e-2 



1 (N+l) - ^aix+y+x-^y- 1 ) 



x 



(34) 
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For odd A, the antisymmetric matrix / has to be sup- 
plemented by an additional row and column with 



oo 

fn,N+i = J dxx'h^e-^+^x™- 1 



2A|( A r +1 )/ 2 _„|(a). 



(35) 



The analytical continuation a — > —ij£e/vF = —ire/N 
can be done directly in the integrals (J54|), ((3^) if accom- 
panied by an appropriate shift of the integration contour 
of the variables x and y in the complex plane. After 
that, numerical evaluation of the partition function j3^ ) 
for complex a and its derivatives, and hence of the DOS 
i^(e), is straightforwardly The small-a asymptote of the 
partition function Z(a) is given by Eq. (^5|) for odd A 
and by 



Z(a) = a~ {N/2)2 [1 + ca In a + 0(a)] 



(36) 



for even N, where c = n[(N/2 - l)!!/(A/2 - 2)!!] 2 if N = 
4n and c = [(A/2 - l)!!/(A/2 - 2)!!] 2 if AT = 4n - 2, 
n = 1,2,. . .. A s a result, for small energies, the DOS has 
the form ( [26a|) for odd A, while 



vqc. 



v{e) = — |ln(£-r)|, 0<er<l, 



(37) 



for even A, where the coefficient c is given below Eq. 
(|36|). In the limit A 3> 1, this asymptote simplifies to 
z/(e) = z/ 1 ln(er)|, < er < 1. 

The DOS is plotted for the three chiral classes and for 
various values of A in Fig. ||. The parity effect and the 
dependence on symmetry class is clearly seen. 



C. BdG universality classes 



We now apply Eqs. (|17|)-(|20|) to the BdG universality 
classes. We choose A = cosh2a;, C = (1, oo), and 



w(X) = exp(-aA)(A 2 - l)^ 1 )/ 2 , 



(38) 



to arrive at the standard form (|l6|). In all calculations 
below we choose the polynomials p„(A) = A" -1 to com- 
pute Z(a). 

When both TR and SR symmetries are present (class 
CI), say as is the case with a dirty ci-wave superconduc- 
tor, we use Eq. (|l^) and find 



Z(a) oc Det 



da 



m+n — 2 



N 



(39) 



With the help of the asymptotic limit a — ■> of the mod- 
ified Bessel function Ki, we obtain 



Z(a) oc a - N(N+r > [1 + ±A(A + l)a 2 In a + 0{a 2 )] 



(40) 



10- 
8 

■£ 6 ' 
> 4- 

2 

0- 




0.005 



0.01 



0.015 





£ T 

FIG. 2. Density of states for the chiral classes. The energy 
is measured in units of 1/r = vf /(N'fl), the DOS is measured 
in units of vq — NcI/tvvf- Figure En, is for the case of preserved 
time-reversal symmetry and spin-rotation invariance (Hamil- 
tonian of class BDI in Cartan's classification), Fig. En is for 
the case of broken time-reversal symmetry (class AIII), and 
Fig. ^c is for the case of preserved time-reversal symmetry 
and broken spin-rotation invariance (class CII). The density 
of states is shown for iV = 1, 2, 3, and 4; the values of N are 
indicated in the figure. 
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Substituting a — > —ire/N and using Eq. ( [14| ) one verifies 
that the DOS vanishes linearly with energy as e — » 

v{e) = (1 + jr) \er\, < er « 1. (41) 

When SR symmetry is present but not TR (class C), 
we find 



Z(a) oc ^ Det [/, 



2N 

™Jm,n=l' 



/mn = (n - m) [ - — 



a \ 1 + a 



(42a) 
(42b) 



In the limit a — > 0, we obtain 

Z(a) oc a-^ 2 ^ [1 + Cl a 2 + c 2 a 3 + 0(a 4 )] , (43) 

where c 2 — j^N(N + X)(2N + 1) and c% is a numerical 
coefficient which is not important for the small energy 
asymptote of the DOS. From Eq. one verifies that 

the DOS vanishes quadratically with energy as e — > 0, 

Ke)^o(l + ^)(l+27v)(£r) 2 , 0< £ r«l. (44) 

The linear and quadratic energy dependencies of the 
DOS in the thermodynamic limit and for small energies 
in the classes CI and C, Eespectively, have been found 
previously by Senthil et a/.tafia The same suppression of 
the DOS appears in the random-matrix theory for finite- 
size systems of classes CI and Cj3 provided the size of the 
system is set equal to a localization volume (correspond- 
ing to a segment of length Nt of the quantum wire). 

When TR symmetry is present but not SR (class Dili) , 
we find from Eq. (0) 



Z(a) oc Det 



d_ 

da 



m+n— 2 



K (a) 



N 



(45) 



Note that N is even for the class DHL In the limit a — > 0, 
using the asymptotic expansion of the modified Bessel 
function Kq, we obtain Z{a) oc g- N( - N -^\lna + 0(1)], 
i.e., corresponding to a DOS of the form (26a). For class 



Dili, the case N = 2 is special: The Lie algebra so(4, C) 
generated by the transfer matrices of a quantum wire in 
class Dili with N = 2 (see Ref. [50|) is not irreducible, 
so(4,C) w sZ(2,C) x sZ(2,C) (see Ref. |). As a result, in 
the case Dili with N — 2 the scaling of the parameters 
Xj and of the eigenphases (f>j is described by two mean 
free paths that need not be equal; one for each copy of 
the non-compact Lie algebra s/(2,M). The scaling equa- 
tions (||) , (g) , and (||) , in which only one mean free path 
appears, correspond to the special case where these two 
mean free paths are equal, which requires fine tuning 
of the disorder. Only in that case the DOS singularity 
In Refs. 



(26a) is found. In Refs. |34] and |35| it is shown that the 
general case of different mean free paths has a power law 
dependence of the DOS for e — > with a non-universal 
exponent. 
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— i 1 — 

0.4 0.6 



0.8 




£ I 

FIG. 3. Density of states for the BdG universality classes. 
Fig. ^a contains results for classes Dill (solid curve) and 
D (dotted curve), for which spin-rotation invariance is bro- 
ken; Fig. is for classes C (solid) and CI (dotted), where 
spin-rotation invariance is preserved. The channel number iV 
is indicated in the figure. 
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Finally, when neither SR nor TR symmetries are 
present (class D), it is enough to use only Eq. ( |l8|) be- 
cause the number of eigenvalues N is a multiple of four 
in this case. We then obtain 



3 



1, . . . , N, where N is the number of independent 



Z(a) oc y/ Det [f m n]*,n=v 



(46) 



fmn 



dx I dy 



x n-l y m-l 



y/(x*-l)(y*-l) 



.m-1 ll-l 



With the help of the asymptotic limit a — > we find 
Eq. ( [26a|) for the asymptotic DOS. Numerical evalua- 
tion of the DOS is possible using @). The substitu- 
tion a — > —ire/N can be done before the integration and 
the convergency of the integrals has to be provided by 
the appropriate shift of the integration contour of x and 
y in the complex plane. A divergence of the DOS was 
also predicted by SenthiLand Fisher for two-dimensional 
systems of class D/DIII.E3 

The case N = 2, which corresponds to a "spin- 
less" class D and was studied in Refs. |34| and [35|, is 
again very special [see the discussion of class Dili with 
N = 2 below Eq. @]: The Lie algebra so(2,2) gen- 
erated by the transfer matrices of a quantum wire in 
class D with N = 2 (see Ref. is not irreducible, 
so(2, 2) w sZ(2, R) x sZ(2, K) (see Ref. g). Therefore, as is 
explained in Refs. p4| and 35, instead of the Dyson singu- 



larity ( 26a ) an algebraic e-dependence with non- universal 
exponent is obtained for the DOS singularity at e — 
whenever the two independent mean free paths associ- 
ated to the two copies of the non-compact Lie Algebra 
s/(2,R) are unequal. 

The DOS is shown for the four BdG classes and for var- 
ious N in Fig. ||[ Note that there is no parity effect for 
the BdG classes and the DOS becomes almost indepen- 
dent of N beyond ./V m 4, provided energy is measured 
in units of 1/r. 



IV. CONCLUSION 

In this paper, we have presented a unified picture for 
both transport and the density of states (DOS) in a 
long quantum wire, for all ten (pure) Cartan symmetry 
classesO While it was well known that the dependence 
of the conductance on the length L of the quantum wire 
could be inferred from the theory of diffusion on symmet- 
ric spaces through the Dorokhov-Mello-E^ejoa-Kumar 
(DMPK) equation and its generalizations ,tjE3E3 we have 
shown that the same geometrical framework underlies the 
dependence of the DOS on the energy e (measured with 
respect to the band center or Fermi energy e = 0). Our 
unified picture consists of a general scheme that permits 
the construction of scaling equations for the eigenvalues 
e 2i0j Q f -j- ne maTTrx product r^(— e)r(e) of reflection ma- 
trices r at different energies for a disordered quantum 
wire and for all ten pure Cartan symmetry classes. [Here 



eigenvalues of (— e)r(e).) Using the analytical contin- 
uation e — > iui, this scaling equation could be obtained 
from that for the reflection eigenvalues Rj of the quantum 
wire in the presence of absorption, which, in turn, can be 
obtained from the DMPK equation. Using the "node- 
counting theorem" that relatea_thje L-dependence of the 
eigenphases cf>j and the DOS,e5l3 we were thus able to 
find the DOS in the thermodynamic limit L — > oo from 
the stationary solution of the scaling equation. 

Relatively little attention has been paid to the behav- 
ior of the density of states (DOS) in the standard prob- 
lem of Anderson localization for a system of infinite size, 
since this DOS, with the exception of its tails, is rather 
insensitive to disorder. This is not true in the presence of 
symmetries for which one energy in the spectrum plays a 
special role, as is the case for the chiral and Bogoliubov- 
de-Gennes (BdG) universality classes, where a sublattice 
symmetry or the presence of a superconducting order pa- 
rameter singles out the band center or the Fermi energy 
e = 0, respectively. In those cases, the DOS is singular 
around zero energy, the energy scale for the singularity 
being h/r, where, for the case of a long quantum wire we 
consider in this paper, r ~ N 2 £/vf is the time for dif- 
fusion through a segment of the wire of length N£. The 
precise form of the singularity depends on the symmetry 
class and, for the chiral classes, on the parity of N. The 
results are summarized in Table || for the limit of a mul- 
tichannel quantum wire, N 3> I. It is in this limit, that a 
universal dependence of the DOS on the symmetry class 
is expected. 

One might be tempted to argue that, for those classes 
where the quasiparticle states are localized at and near 
zero energy, the singular behavior of the DOS in the ther- 
modynamic limit is controlled by the (zero-dimensional) 
random matrix theory for an effective finite-size system 
with the linear size of the localization length £ = Nl. 
Whereas this argument predicts the correct energy scale 
H/t for the singularity, it fails to reproduce the extra log- 
arithm in the functional form of the singularity for the 
chiral classes, see Table 0. The origin of the logarithm 
in the DOS can be found in the presence of two com- 
peting localized modes in that case, that each have the 
same localization length, see Refs. |2^^l|. As was shown 
in Ref. [H], the presence of a small amount of dimeriza- 
tion of the hopping amplitudes lifts the degeneracy of 
the localization lengths and removes the logarithm from 
the functional form of the DOS singularity. On the other 
hand, for the B4G classes C and CI, there is only one 
localized modej^j which explains the absence of the ex- 
tra logarithm in the functional form of the DOS in those 
classes. 
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class 


TR 


SR 


v(e) for < er < 1 




Yes 


Yes 


!A)|ln(er) 


chiral, TV even 


No 




7r^ |c:Tln(er)| 




Yes 


No 


(7r^o/3)|(£r) 3 ln(er)| 


criirfil /V nrlH 

V -* 1 ,1 1, 1 CL1 -. 1 V W v I \_L 






li l/y / C / 111 1 C / 1 




Yes 


Yes 


(nvo/2)\er\ 


BdG 


No 


Yes 


vo\et\ 2 




Yes 


No 


■kvq/\et ln 3 (eT)| 




No 


No 


■kvqI\et ln 3 (£-r)| 



TABLE II. Asymptotes of the DOS v(e) for small energies 
< er 1 and for N S> 1, where r = N-jI/vf, 7 is given in 
Eqs. © and ©, £ is the mean free path, and i>f the Fermi 
velocity. (Expressions for the DOS asymptotes for finite N 
can be found in the text.) 
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